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Abstract 

Within the lowest-order Born approximation, we present an exact calculation of the time dy- 
namics of the spin-boson model in the Ohmic regime. We observe non-Markovian effects at zero 
temperature that scale with the system-bath coupling strength and cause qualitative changes in 
the evolution of coherence at intermediate times of order of the oscillation period. These changes 
could significantly affect the performance of these systems as qubits. In the biased case, we find a 
prompt loss of coherence at these intermediate times, whose decay rate is set by ^/a, where a is 
the coupling strength to the environment. These calculations indicate precision experimental tests 
that could confirm or refute the validity of the spin-boson model in a variety of systems. 
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Novel solid state devices that can control spin degrees of freedom of individual 
ele*o,4 a. o. d.c.e.e .a.e. in .pe.conducting ccit^B Q Q Q, ^How 

promise in realizing the ideal of the completely controllable two-state quantum system, 
weakly coupled to its environment, that is the essential starting point for qubit operation in 
quantum computation. From a fundamental point of view, these experimental successes also 
bring us close to embodying the ideal test of quantum coherence as envisioned by Leggett 



many years ago 



7[, in which a simple quantum system is placed in a known initial state, is 
allowed to evolve for a definite time t under the action of its own Hamiltonian and under 
the influence of decoherence from the environment, and is then measured. 

Rece,. e.pe.,„e„ts. starting w.th Q. sKow that tKi. idea, test can be ..npie.ented in 
practice. The decay of quantum oscillations due to environmental decoherence is nowP, 0, 
la |6[ sufficiently weak that some tens of coherent oscillations can be observed. If quantum 
computation is to become a reality, it is believed ^ that these systems will eventually need 
to achieve even lower levels of decoherence, such that thousands or tens of thousands of 
coherent oscillations could be observed. This prospect of producing experiments with ultra- 
long coherence times in quantum two state systems offers a new challenge for theoretical 
modelling of decoherence. Despite the many years of workj^,^^ following on Leggett's initial 
proposals, there has never been a full, systematic analysis of the most popular description of 
these systems, the spin-boson model, in the limit of very weak coupling to the environment. 

In this Letter we provide an exact analysis of the weak coupling limit of the spin boson 
model for the Ohmic heat bath, and in the low temperature limit. In this limit the Born 
approximation (to the self energy) should become essentially exact, and we make no other 
approximations in our solutions — in par ticular, no Markov approximation is made. As 
other workers have recently emphasized 3, Q| , understanding the details of the short-time 
dynamics of this model is especially crucial for the operation of these systems as qubits. 

We find important, new, non-Mar kovian effects in this regime. At lowest order in the 
Born expansion of the self energy superoperator, the time dynamics of the model rigorously 
separates into a sum of strictly exponential pieces (the usual "Ti" and "T2" decays of the 
Bloch-Redfield model) plus two distinct non-exponential pieces that arise, technically speak- 
ing, from two different kinds of branch cuts in the Laplace transform of the solution of the 
generalized master equation that we obtain. 

These two contributions both have power-law forms at long times, t > Ti,T2, and thus 
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formally dominate the exponentially-decaying parts. But more interesting is that they both 
give new structure to the time evolution at intermediate times t, l/uc < t < Ti,T2] this 
structure typically occurs for t on the order of the oscillation period. (Here, Uc is a high 
frequency cut-off of the bath modes, defining the very short time regime, t < l/oOc, which is 
of no interest here.) We can explain our results in the language of the double-well potential, 
where the two quantum states are "left" and "right" (L/R), the t = state is pure L, and 
the system oscillates in time via tunneling from L to R. The first branch-cut contribution 
is most important in the unbiased case (L and R energies degenerate) and it causes the 
system, starting immediately in the first quantum oscillation, to spend more time in the R 
well, that is, the opposite well from the one the system is in initially. The second branch- 
cut contribution, present when the system is biased, adds to the amplitude of the coherent 
oscillation, but dies out after an intermediate time which scales like the inverse square root 
of the interaction strength a with the bath. This prompt loss of coherence, whose amplitude 
is proportional to a, changes qualitatively the picture of the initial decay of coherence that 
is so important for discussions of fault-tolerant quantum computation . 

We are interested in studying the time dependence of the system density matrix ps{t) = 
TTBp{t) with a time-independent system Hamiltonian, and in the presence of a fixed coupling 
to an environment. An exact equation for ps - the generalized master equation (GME) - 
isQ 

Ps{t) = -iLsPsit) - I C dt'J^sit - t')ps{t'), (1) 

Jo 

J^sit) = -tTTBLsBe~"^'''LsBPB- (2) 

Here the kernel ^s{t) is the self energy superoperator, the system-bath Hamiltonian is 
written H = Hs + Hsb + Hb, the Liouvillian superoperator is defined by L^p = [Hx,p], 
Pb = e^^^B /Z, 13 = l/fc^T, T is the temperature, and Q is the projection superoperator 
Q = 1 — pbTtb- Eq. ((T)) is written for the case TtbHsbPb = 0, and the total initial state 
is taken to be of the form p(0) = ps{0) ® Pb- Since we are interested in the case of weak 
coupling to the bath, we will consider a systematic expansion in powers of this coupling Lsb 
m the self-energy operator E.(t). Rete,rt,o,r of only the lowest order term .n this expansion, 
giving the Born approximation, is obtainedjla] by the replacement e e '^Qi^s+LB)t 

Eq. ©. 

We now proceed to solve the GME with no further approximations. This distinguishes 
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our work from previous efforts, in which various other approximations 



LS fsecular, rotatin; 



wave, Markov, "non- interacting bhps", short time) are made (see, e.g., |7l. l9l llOL IllL Il2l|). 
We will find that, in particular, avoidance of the Markov approximation endows the solution 
with qualitatively new features. 

We obtain our solution for the special case of a two-dimensional system Hilbert space, 
and a system-bath coupling of a simple bilinear form, Hsb = S^X {S (X) is an operator in 
the system (bath) space). In this case the GME in Born approximation can be rewritten 
in an ordinary operator form: 

{&^{t)) = -tTTsa,[Hs, ps{t)] - f dt'I^it, t') , (3) 

<j 

i,{t,t') = i^,{t') + j2i,At')Mt~t')), (4) 

u=l 

V(t') = Re {Ci-t')TTsa,i-t')[a^, 5]5(-t')} • (5) 

Here (Jx,y,z are the Pauli operators, ctq = I, {x) = Tr sxps, the bath correlation function is 
C{t) = TTB[XX(t)pB] = C'(t) +iC"(t), the time dependent operators are in the interaction 
picture, i.e., o{t) = Q'^i^s+HB)tQ^-i{Hs+HB)t ^ ^^^^ (j/ g^^^^ (jn fjgnote the real and imaginary 
parts of the bath correlator. 

Without loss of generality, we can take the system operators to be of the form 

A e 

Hs = -^(Tx + ^(Tz (6) 

and S = az- Then the GME can be written in an explicit form {&.{t)) = R* {a} + k. Here 
a denotes the vector {ax, cr^, a^)"^, convolution is denoted A* B = /q dt' A{t')B[t — t'), and 

^ -firi(t) -eS{t) + %K^{t) ^ 
R(t) = e6{t) - lK;{t) -Tyit) -A5{t) , (7) 
^ A(5(t) J 

m = (-|A;-(t), -k;{t), (8) 

with E = + A2, and kernels given by ri(t) = ^ cos{Et)C'{t), Ty{t) = ^(1 + 
^cos{Et))C'{t), K+{t) = ^sin(Et)C"(t), k-{t) = ^ dt' sm{Et')C"{f), and k;{t) = 
^ /o* dt'il - cos{Et'))C"{t')\i. 

These equations can be solved in the Laplace domain. Defining the Laplace transform as 
f{s) = e~'^^f{t)dt, the solutions are, for the "standard" initial conditions {g_{t = 0)) = 



(0, 0, 2;o = 1 



/ / N\ , , AN(s) zo 

N{s) = I - ^K(')) ^"(^) + (7^0 + Kis)) (3 + |^ri(.)^ , (10) 
D{s) = + r,(.) + v) + S^^^')) + " f ^^""('O' • ^^^^ 

We can give more explicit solutions if we specialize to the spin-boson model, for which 
^ = I]nCn(^Ji + bn) and Hb = J2n^nbnbn- If the spectral function J^w) is defined as 
J (u) = J2n els {u -ujn), then 

C{t)= duJ{u){coth{l3u/2)cos{ut) + ism{ujt)). (12) 

JO 

For most of the sequel we will consider the Ohmic case, for which J{uj) = ^cue^^^'^'' {cUc is 
an ultraviolet cutoff). In this case Eq. (fT^ becomes 



where ifj' is the derivative of the digamma function|l7l|. 

For discussing the exact solution it is instructive to understand the structure of the 
solution in a Markov approximation. This approximation is obtained by replacing all the 
kernels Fi, F^, Ky, k^, and ky by their values at s = 0. In this case the solutions Eqs. 
Q are rational functions of s. Then, if the poles of these rational functions are located at 
positions in the complex s plane, with residues rk/27ri, then the inverse Laplace transform 
can be written (cr^(t)) = J2k^k ^wi^it). We indicate here that while the residues do depend 
on the label /i = x, y,j, the pole positions do not, as is suggested by the form of Eqs. Q. 

As is well known[9,], there are four poles: Si = 0, S2 = — F?, and 53,4 = — F^ ± iE. 
The first pole describes the long-time asymptote of the solution (stationary state), the 
second the purely exponential, "Ti"-type decay (relaxation), and the last two (complex 
conjugate paired) describe an exponentially decaying sinusoidal part, the "T2"-type decay 
of coherent oscillations. The residues of these poles are, to lowest order in a, = x^o = 
-{A/E) tanh(/3E/2), rf = = 0, = = -{e/E) tanh(/3E/2), = eA/E^ - x^, 
r| = 0, r| = eV^^ - z^o, r^4 = -tA/2E'^, = -A/2E, and rl^ = A'^/2E'^. The pole 



positions are, again to lowest order in a, given by = ^ = coth(/?£'/2), = T2 ^ = 
ry2 + ^ffkeT, and E = E + 6E, 6E = dE^^""^ + 6E^'^'^ll8], SE^^""^ = s|i(C + ln(E/cUc)) 
and SE^^^"^^ = ^^(Re ip{iEP /2'n-) — ln(£'/9/27r)), where we have dropped terms of order 
E/uJc and higher, C is the Euler constant, and is the digamma function These 
expressions are straightforwardly derivable, and agree with the literature^, except for the 
energy shift due to vacuum fluctuations, ,5£;Lamb (^-^j-^j^j-^ contains in general ln{E/uJc) and not 
ln(A/co'c)). We note that this Markovian theory satisfies the expected fundamental relation 
FO = FO/2 + {2eyE^) J^^dt{X{t)X)B (Korringa relation)0; 

also, to lowest order in a, 

the asymptotic values of (cr^(t 00)) go to the Boltzmann equilibrium distribution of the 
system, e.g., z^o = —(e/E) tanh(/3i?/2), unlike, for example, the popular "non-interacting 
blip" approximation [9]. 

We now return to the exact solution, examining it in detail at vanishing temperature 
T = 0. In this case the Laplace transform of C in Eq. ()13|) is 

Ct=o{s) = q;s/2 (— cos (s) Ci (s) — sin (5) si (5)) 

— ia/2 {—Uc + s sin (5) Ci (5) — s cos (5) si (5)) , (14) 



where s = s/uJcll9\. There is an important feature of this correlation function that makes the 



Markov solution qualitatively incomplete: whi 



e the sine integral si is analytic, the cosine 



integral Ci(s) behaves like ln(s) for s 17[. This means that C{s) is nonanalytic at 
s = — it has a branch point there. Thus, the exact solutions (cr^(s)) have extra analytic 
structure not present in the Markov approximation, and the real-time dynamics has 
qualitatively different features in addition to the pole contributions we have just discussed. 

The s = branch point in C{s) leads the kernels Fi(s), Ky{s), and k~{s) to have branch 
points at s = ±iE; the kernels Ty{s) and ky{s) have three branch points, at s = and 
s = ±iE. Thus, the solutions to the GME {(Tx^y^zis)) also have three branch points in the 
complex plane. We find by numerical study that the exact solutions still have four poles 
as before, which, for small a, have nearly (but not exactly) the same pole positions and 
residues as in the Markov approximation. 

Thus, the structure of the solutions in the complex s plane is as shown in Fig. The 
locations of the branch cuts are chosen for computational convenience, as discussed shortly. 
Given this branch-cut structure, the inverse Laplace transform (the Bromwich integral) is 
evaluated by closing the contour as shown. Thus, the exact inverse Laplace transform can 



be expressed as (t > 0) 

Mt)) = -L / dse^\a,is)) = -i. / dse^\a,is)) 
Zm jc Zm JCo 

1 ^ 

-7r-J2^k dxe'i''''\{a^{qkX + r]k)) - {(y^{qkX-r]k))). (15) 

Here = e*^* and r]^ = r]e^^^''~'^^'^\ with i] an infinitesimal positive real number. That is, 
?7fc is an infinitesimal displacement perpendicular to the direction of branch cut k. For the 
cut choices we have made, 9i = 57r/4, 02 = tt/2, 6^ = 37r/2, pi = 0, and P2 = Ps = E. The 
closed-contour integral in the expression can be written as a sum over the four poles, and so 
gives complex exponential contributions to the solution as in the Markovian case. The extra 
terms, the sum over the three branch cuts, are new and give qualitatively different features. 
The contributions of the second and third branch cuts are complex conjugates of each other, 
so we will discuss them together as the "branch cut 2" ("bc2" for short) contribution. The 
first cut will be discussed as the "bcl" contribution. 

The contribution of these cuts to the solution is independent of the detailed positioning 
of the branch cuts, so long as they are not moved across a pole; the choice of the direction 
of bcl is a computational convenience — the apparently most natural choice of this cut 
direction, along the negative real axis, passes it essentially on top of the Fi pole, making 
the evaluation of the branch-cut integral numerically inconvenient. As a check, we find that 
the results we discuss now are indeed independent of the cut direction. 

We have done a detailed study of these branch-cut contributions for {az{t)) = z{t). For 
the unbiased spin-boson case, e = 0, an essentially analytic calculation can be done for all 
contributions. In this case there is no bc2 contribution, Zi,c2 = 0. The bcl contribution can 
be obtained analytically to leading order in a: z(t) = Zpoiesit) + -2bci(^), 

z^^^{t) = -a{l - At[Ci(At) sin(At) - si(At) cos(At)]}. (16) 

This function, plotted along with the pole contribution in Fig. for the choice of parameters 
shown, has the following features: 2feci(^) is negative for all t, it is monotonically increasing, 
and its long-time behavior is Zbciit) ~ — 2a/(At)^. Also, Zbciit = 0) = —a. 

Let us survey, then, the peculiar features that this branch cut contribution introduces 
into the time response z{t). Visualizing the e = spin-boson model as a symmetric double 
well system coupled to its environment, the bcl piece being negative means that, if the 
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system is initially in the left well, it will, in the course of coherently tunnelling back and 
forth, spend more time in the right well! This effect becomes strongest at long time, much 
longer than T2, for in this regime the pole contributions are exponentially small, while the 
bcl contribution decays like a power law. Experimentally it may be hard to see the effect in 
this regime (on account of finite-temperature effects, for example), so it is important to note 
that this memory effect appears already at early times, indicating that already in the first 
couple of coherent oscillations, there will be an excess amplitude in the right-well excursions 
as compared with the left- well excursions, by an amount proportional to a. We believe, 
on the basis of a variety of evidence 2l|] , that the Born approximation should be reliable 
up to a's of order 1 — 2%; thus, experiments that look at coherent oscillations accurately 
at the percent level (which, it seems, will ultimately be necessary for performing quantum 
computation) could readily see this bcl effect. 

We note several other interesting features of our solution for e = 0. Taking into account 
the non-Markovian effects, we can do a more precise calculation of the pole positions and 
residues (only poles 3 and 4 contribute). We find, for T = 0, r2 = — Re(s3) = F^r, where, 
as before = air A/2, and the renormalization factor r is given by r = (1 — a)/{H^ + 
aV^) < 1, with K = I- 2a(l/2 + C + ln(A/cjJ). Further, Im(s3) = E + SE^^'^^f, with 
f = (k — a7r^/2(C + ln(A/co'c)))/(K^ + (avr)^). These expressions are obtained as systematic 
expansions in the small parameters T2/E and 6E/E, and they match a direct numerical 
evaluation of the pole positions very well up to a's of a few percent. For the corresponding 
pole residues we find the simple result in leading order + = 1 + a + 0(q;^). This would 
be impossible in a Markovian theory, in which z{t = 0) = + r^, so that + would 
be exactly 1 to all orders in a. In fact this excess pole residue is exactly what is needed to 
cancel out the initial value of the bcl contribution to z{t). We note that our results for the 
residues differ from the weak-coupling expressions in the literature (we are not aware of 
prior reports on r, f). 

For the biased model (e 7^ 0) the bc2 contributions become nonzero; we find that they 
give other peculiar non-exponential corrections to the solution z{t), very different from the 
bcl contribution. We can do a nearly analytic evaluation of the bc2 contribution to Eq. 
()15|): Using Eq. Q and expanding to lowest order in a, we find for the integrand of the 
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sum of the k = 2 and 3 terms of (fT3j) . 

LJ [E^ — LU^ + b+{uj))^ + [uj)^ 
Here b{iuj ± r]) = b'^{uj) ± ib^{uj), b{s) = a{d{s) + n{s){s'^ + E'^)/A), where d{s) and n{s) 
are given by N{s) = A + an{s) (see Eq. (UHl)) and L'(s) = + E'^ + ad{s) (see Eq. (ITT|) 'l. 
Since b~{uj) = for |ti;| < E, it is reasonable to expect that b~ will grow linearly as one 
passes onto the branch cut; and, in fact, we find from numerical study that a good ansatz 
is b^{uj) = {E ~ uj)b^{uj), with b^{uj) being a weakly varying, real function of u/E. With 
this, for uj of order E, Eq. (|T7jl simplifies to 

. A%~iE) 1 

Zbc2is = lUj) ~ ^ . (18) 

^ ' 2E^ u~E ^ ^ 

We find that (fTH|) should be valid for uj > E + b^{E) /2E. Using (fT5|l we can do the branch 
cut integral, which gives (for t < l/{aExQ) — see Appendix for an alternative approach), 

Zbc2{t) ~ <yxi \og{xoaEt) cos{Et + 0). (19) 

Here the dimensionless constants xq = \b'^{E)\/2aE'^ = \6E\/aE and xi = A%~'{E)/2aE^. 
Since 6^ oc a, these constants are independent of a. The last expression for xq comes from 
an evaluation of b^{E): it is directly related to the energy renormalization in the Markov 
approximation, &+(£') = 2ESE^'''^^. 

In Fig. 121 we show a direct numerical evaluation of Zi,c2{t). One can see the decay of the 
oscillatory part, which is logarithmic according to Eq. (fT^ . Even though the decay is very 
non-exponential, it is reasonable to attempt to characterize this decay by a time scale. Eq. 
()19p obviously does not work at t = 0, since it is logarithmically divergent. This is not 
surprising, since our calculation has neglected cutoff effects (dependence on cj^), so Eq. (fTI?|) 
is not expected to be correct for t < l/oOc- However, if we consider "early" time to be the 
first half-period of the coherent oscillation, to = ^^/E, then Eq. ()19|) should be valid and we 
can use it to characterize the decay by determining the time th at which Zbc2{t) decreases to 
half its early-time value, i.e., Zbc2(th) = hzbc2(to)- We obtain 



1 vrE 1 1 
''=EiW\"'E7^- 
Surprisingly, th oc l/y/a depends non- analytically on a. This explains the effect that is 
evident in Fig. |21 for small a, th <^ T2, that is, on the scale of T2, there is a very rapid 



loss of coherence as contributed by bc2. This phenomenon may be called a prompt loss of 
coherence, as it would appear experimentally as a fast initial loss of coherence (from 100% 
to (1 — ca)100%, c being some constant near unity), followed by a much slower, exponential 
decay of coherence on the regular T2 time scale. 

We give a few final remarks about the bc2 calculation. The absolute size of the bc2 
contribution reaches a maximum near the value of e/A used in Fig. |21 the relative size of 
this contribution continues to increase as |e|/A increases, so that it eventually becomes much 
larger than the pole contribution (but all contributions to z{t) go to zero as |e|/A — 00). 
When |e| ~ A, we find that, because of the prompt loss of coherence, there is a deficit in the 
total pole contribution, that is, J2k''^k = 1 " 0{a) < 1. Even in the absence of an explicit 
branch cut computation, this deficit signals the prompt loss of coherence, in that it indicates 
that the exponentially decaying contributions to z{t) do no account for all the coherence 
near t = 0. Note that this is opposite to the unbiased case, where, as a result of the bcl 
part, there is an excess pole contribution. 

Naturally, many more regimes could be studied using the present approach. Recently, 
there has been interest in varying both the svstemjl^ and bathjs^ initial conditions, as 
well as in varying 

the systematic Born expansion procedure we report here can be done. It is clear on general 
grounds that the appearance of branch cut contributions will not be restricted to the Ohmic 
model we have studied in detail here; It is easy to show that for any spectral density of the 
form J{uj) oc cu" at low frequencies (n = 0, 1, 2, ...), C(t) will have a power-law dependence at 
long time, and thus C{s) will have a branch point at s = 0. So, interesting non-exponential 
contributions to the dynamics are expected in all these cases. 

Our hope is that, using the present and further exact calculations of the weak-coupling 
behavior of the spin-boson model, a tool will be made available to permit precision ex- 
periments to test the validity of the model (which, at present, is only phenomenologically 
justified) in various physical situations of present interest in quantum information. A fun- 
damentally correct, experimentally verified theory of the system and its environment should 
ultimately be of great value in finding a satisfactory qubit for the construction of a quantum 
information processor. 

DPDV is supported in part by the National Security Agency and the Advanced Research 
and Development Activity through Army Research Office contract number DAAD19-01-C- 



10 



0056. He thanks the Institute for Quantum Information at Cal Tech (supported by the 
National Science Foundation under Grant. No. EIA-0086038) for its hospitality during the 
initial stages of this work. DL thanks the Swiss NSF, NCCR Nanoscience, DARPA and the 
ARO. 



[1] R. Hanson et al, "Zeeman energy and spin relaxation in a one-electron quantum dot," 
|cond-mat/0303139, 

[2] D. Loss and D. P. DiVincenzo, "Quantum Computing with Quantum Dots," Phys. Rev. A 
57, 120 (1998), |cond-mat/9701055 



[3] Yu. Makhlin, G. Schon, and A. Shnirman, "Quantum state engineering with Josephson- 

junction devices," Rev. Mod. Phys. 73, 357-400 (2001), |cCTl d-mat/00 11269| 
[4] Y. Nakamura, Yu. A. Pashkin, and J. S. Tsai, "Coherent control of macroscopic quantum states 

in a single-Cooper-pair box," Nature 398, 786-788 (1999); Yu. A. Pashkin, T. Yamamoto, O. 

Astafiev, Y. Nakamura, D. V. Averin, and J. S. Tsai, "Quantum oscillations in two coupled 

charge qubits," Nature 421, 823-826 (2003). 
[5] I. Chiorescu, Y. Nakamura, C. J. P. M. Harmans, and J. E. Mooij, "Coherent Quantum 

Dynamics of a Superconducting Flux Qubit," Science 299, 1869-1871 (2003). 
[6] D. Vion et al, "Manipulating the quantum state of an electrical circuit," Science 296, 886-889 

(2002). 

[7] A. J. Leggett et al, Rev. Mod. Phys. 59, 1 (1987). 

[8] D. P. DiVincenzo and D. Loss, "Quantum Computers and Quantum Coherence," J. of Mag- 
netism and Magnetic Matls. 200, 202 (1999), cond-mat/9901137 l 
[9] U. Weiss, Quantum Dissipative Systems (World Scientific, Singapore, Second Edition, 1999). 
[10] M. Grifoni and P. Haenggi, "Driven Quantum Tunneling," Phys. Rep. 304, 229-358 (1998). 
[11] Daniel Braun, Fritz Haake, and Walter T. Strunz, Phys. Rev. Lett. 86, 2913 (2001) 
lquant-ph/0006117); Walter T. Strunz, Fritz Haake, and Daniel Braun, Phys. Rev. A 67, 
022101 (2003) ( .quant-ph/0204129 1. 



[12] V. Privman, "Onset of Decoherence in Open Quantum Systems," cond-mat / 0303157} 



"Initial Decoherence of Open Quantum Systems," J. Stat. Phys. 110, 957-970 (2003), 
quant-ph/02 05037| 



11 



[13] J. Preskill, "Fault-tolerant quantum computation," in Introduction to Quantum Computation 
and Information (ed. H.-K. Lo, S. Popescu, and T. Spiller, World Scientific, Singapore, 1998), 
pp. 213-269. 

[14] E. Fick and G. Sauermann, The Quantum Statistics of Dynamic Processes (Springer, Berlin, 
1990). 

[15] M. Celio and D. Loss, Physica A 150, 769 (1989). 

[16] Note that in this theory there is a very simple relationship between ay and cr^: Aay = &z- 
[17] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, (Dover, New York, 
1965), Chaps. 5 and 6. 

[18] C. W. Gardiner and P. Zoller, Quantum Noise (Second Edition, Springer, Berlin, 2000). 
[19] The constant term uJc in C^^q{s) cancels in the physical correlators k^{s) and ky{s). 
[20] A. Abragam, The Principles of Nuclear Magnetism (Clarendon Press, Oxford, 1961), pp. 
2.12-13. 

[21] For example, the scaling relations discussed in the Appendix work very well up to a ~ 0.01. 
[22] A. Shnirman, Yu. Makhlin, and G. Schon, "Noise and decoherence in quantum two-level 
systems," |cond-mat/0202518| 

APPENDIX A: SCALING FORM FOR BRANCH CUT INTEGRALS 

By numerical study we find that the branch cut integrals conform to some simple scaling 
laws for small a. If we write the bcl and bc2 integrals as zi)ci{t) = dxe~'^^^^ zi,ci{s = qix) 
and Zbc2it) = Re dxe''^^ Zt,c2{s = ix), then we find that for small a and for s << Uc, Zbci{s) 
can be written in a scaling form 

zi>,,{x) = {a/E)f,{e/A,x/E). (Al) 

But for bc2 a very different scaling law applies: 

Zbc2{x) = (l/E)/2(e/A, {x - E)/aE). (A2) 

Here /i^2 are dimensionless, "universal" functions that govern the behavior of the branch cut 
contributions for small a. For bcl, the behavior that the scaling law gives is very simple: Eq. 
()A1|) implies that Zhciif) = afi{e/ A, Et), where /i is the Laplace transform of the scaling 
function /i. We might have expected this behavior from Eq. ()16|). from which we can read 
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off tlie scaling function for e = 0. In fact it appears from numerical studies that /i hardly 
changes as e is varied, except for an overall scale factor; that is, /i(e/A, Et) ^ a{e/ A)b{Et). 
We find that the scaling function a(r) > is peaked at r = 0. So, the memory effect 
described above for e = persists for finite e, but becomes smaller. For |e| ^ A the bc2 
contribution, which we will describe now, becomes dominant over the bcl one. 

Returning to Eq. ()A2|) . if we write the Fourier transform of the scaling function as 
/2(r) = /o°° e*^'^/2(x)(ix and consider its polar form /2(t) = r2(r)e*'^^'''^\ then we obtain 



This shows that bc2 contributes an oscillatory part to the solution, whose "T2" decay is 
determined by the features of the scaling function r2. A few more observations about /2 
(obtained initially from numerical study) reveal some crucial properties of the r2 function: 
1) /2(0) = 0; 2) 1/2(3^)1 has a single maximum at x = xq, where xq is some constant of order 
unity; 3) Most important for the present discussion, for x > Xq /2(x) approaches 1/x, that 
is, f2{x) ~ Xi/x, where Xi is another real constant of order unity. Fact 3) implies that, for 
r ^ 0, T2(r) ^ Xi log(xoT). That is, we conclude that at sufficiently short time (actually for 
t < l/{aExo), so a relatively long time). 



Zbc2if) = Oir2{oiEt) cos{Et + 02(tt-Et)). 



(A3) 



Zbc2(t) = axi \og{xoaEt) cos{Et + 0) 



(A4) 



as stated in the text. 
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FIG. 1: (a) Structure of the solutions in the complex s plane. The four poles pi, p2, ps, 

and p4, are indicated by crosses: the three branch points at s = 0, ztiE are indicated by solid 
circles, and the three branch cuts chosen, bcl, bc2, and bc3, are indicated by dashed lines. The 
inverse Laplace transform requires an integration along the contour C parallel to the imaginary 
axis. This integral may be evaluated by closing with a contour in the left half plane (Co, the 
Bromwich contour), which lies at infinity except for looping back around each of the branch cuts, 
(b) Zpoiesit) and zi,ci{t) for the unbiased case, e = 0, A = 1, a;c = 30, T = 0, and a = 0.01. t is in 
units of l/E (i.e., E = 1). 
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FIG. 2: Zpoies{t) and 2:fcc2(*) for the biased case, illustrating the prompt loss of coherence produced 
by bc2. Here E = 1, e/A = —1.38, uJc = 30, T = 0, and a = 0.01. For these parameters, the time 
scale for the prompt loss of coherence (using Eq. 1)201) ) is th = 18.98. th is the time at which the 
envelope of Zbc2 falls to half its value at to = ir/E. This time scale is much shorter than the regular 
exponential decay of coherence in ZpoZes! for our parameters, T2 = 204.6. 
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